Entangled-photon generation in nano-to-bulk crossover regime 
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We have theoretically investigated a generation of entangled photons from biexcitons in a semi- 
conductor film with thickness in nano-to-bulk crossover regime. In contrast to the cases of quantum 
dots and bulk materials, we can highly control the generated state of entangled photons by the 
design of peculiar energy structure of exciton-photon coupled modes in the thickness range between 
nanometers and micrometers. Owing to the enhancement of radiative decay rate of excitons at 
this thickness range, the statistical accuracy of generated photon pairs can be increased beyond the 
trade-off problem with the signal intensity. Implementing an optical cavity structure, the generation 
efficiency can be enhanced with keeping the high statistical accuracy. 

PACS numbers: 42.65.Lm, 42.50.Nn, 71.35.-y, 71.36.+C 



Entangled photons play an important role in quantum 
information technologies, and the pursuit of high-quality 
generation of them has become an active research area 
in the fields of quantum optics and solid-state physics. 
In addition to the standard method of generating the 
entangled-photon pairs by using the parametric down- 
conversion (PDC) in second order nonlinear crystals [ij, 
the generation scheme using a semiconductor quantum 
dot [2|, y] attracts much attention because we can gener- 
ate a pure single pair of entangled photons that is highly 
desired in quantum information processing. On the other 
hand, the development of entangled photons as an exci- 
tation light source is becoming of growing importance for 
the next-generation technologies of fabrication and chem- 
ical reaction For this purpose, extra high-power and 
high-quality entangled-photon beams are absolutely nec- 
essary. However, any schemes that lead to the generation 
of such photon beams have not been found yet, although 
several generation schemes have been investigated so far. 
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FIG. 1: Sketches of entangled-photon generation from biex- 
citons (RHPS). (a) In a bulk semiconductor, statistical accu- 
racy worsens by the nonradiative decay of generated photons, 
(b) In a nano-film, most of photons can go outside by the ex- 
citon superradiance, and better accuracy is expected. Center- 
of-mass motions of excitons and biexcitons are confined in the 
film. 



For the high-power generation of entangled photons, 
one of the promising candidates is the method involv- 
ing the biexciton-resonant hyperparametric scattering 
(RHPS) demonstrated in a CuCl bulk crystal This 
generation method is highly efficient owing to the giant 
two-photon absorption in CuCl crystals 0, @|- However, 
as indicated in the experiments (a, [6| , we must consider 
not only the generation efficiency but also the statistical 
accuracy, i.e., accidentally emitted photon pairs (noise), 
which come from independent biexcitons and have no 
entanglement (see Fig. 1(a)). In the scheme for high- 
power generation of the entangled photons, the amount of 
entangled pairs (signal) and signal-to-noise (S/N) ratio 
are always a trade-off, because the nonradiative damp- 
ing of excitons significantly worsens the S/N ratio of 
the entangled-photon generation. Although this trade- 
off problem is seemingly inevitable by the exsistence of 
the nonradiative damping under the resonant excitation 
of the elementary excitation in condensed matters, we 
have revealed that it can be overcome by the enhance- 
ment of radiative decay rate of excitons in the crossover 
regime between the quantum dots and bulk materials . 
Further, we can also control the characteristics of gener- 
ated photons by a variety of degrees of freedom in con- 
densed matters and by the peculiar energy structure in 
the crossover regime 
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Our calculation method is based on theguantum elec- 
trodynamics (QED) theory for excitons 12|. Compared 
to the previous theories |13| , we have developed a frame- 
work applicable to nanometer and submicrometer films 
including explicit degrees of freedom of the center-of- 
mass motion of confined excitons. We consider excitons 
as pure bosons and introduce an exciton-exciton interac- 
tion into the system to discuss the RHPS process. The 
Hamiltonian of the excitonic system is represented as 
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V, 



(1) 



where is the annihilation operator of an exciton in 
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state /i, Up is its eigenfrequency, and 



is the 



nonlinear coefficient. The equation of excitons' motion 
is derived in uj (frequency) domain as f 4 | 
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Here, E + {r,uj) is the positive-frequency Fourier trans- 
form of the electric field at position r. "P^(r) is the coef- 
ficient of excitonic polarization P ex (r) = E M "P ' ^{r)b^ + 
H.c., and is represented as 7 ? J u(r) = Mgr^(r) with the 
transition dipole moment M and the exciton center-of- 
mass wavefunction <?^(r). £> M (u;) is the fluctuation oper- 
ator due to the excitonic nonradiative damping with rate 
7 0X . The equation of motion of E + (r, u) is derived as 

V x V x E+(r,uj) - (uj/c) 2 e hg (r,u)E + (r,u) 
= ifi ujJo{r,uj) + fi u 2 P^(r,uj), (3) 

where Shg(r^uj) is the background dielectric function, /xo 
is the magnetic permeability in vacuum, and Jp(r, to) 
governs the fluctuation of the electromagnetic fields [12| ■ 
The third term on the right-hand side of Eq. ^ is the 
nonlinear term for the entangled-photon generation, and 
it is approximately treated as follows. First, we calcu- 
late the amplitude = (b^(uj)) of excitons in linear 
regime from Eqs. @ and (0) by neglecting the nonlinear 
term and by assuming a pump field as a homogeneous 
solution of Eq. in classical framework of light. Then, 
the biexciton amplitude B\{lj) is calculated from /3 M (w) 
and a phenomenologically introduced damping rate 7b x 



of biexcitons LJ] . Considering a sufficiently strong pump 



power compared to the vacuum fluctuation, we can ap- 
proximately rewrite the nonlinear term in Eq. ^ as 
Ea,,(^ + t^u ~ hQ x )F Kll , v da/ bt(u ~ lo')B x {lo>), 
where fix is the eigen frequency of biexcitons in state 
A and Fx ^ v is the coefficient giving a biexciton state 
l A > = \ E^^W&J&JIO). While fix and Fx <fl<v should 
be determined from Eq. ([T]), we instead assume them 
as follows. If the center-of-mass motion of the low- 
est biexciton level (zero angular momentum) is con- 
fined in a crystal sufficiently larger than its Bohr ra- 
dius, we can approximately express the coefficient as 
Fx,fi,u = Sx.^ Jdr G x {r)g*{r)gt{r), where \<P\ 2 rep- 
resents the effective volume of relative motion of biexci- 
ton, 5x,[i,i> describes the polarization selection rule in the 
excitation and collapse of biexcitons, and Gx(r) is the 
center-of-mass wavefunction of biexciton state A. 

Under the above approximation, by simultaneously 
solving Eqs. @ and ([3]) up to the lowest order, we 
evaluate correlation functions of E + (r,u) from com- 
mutation relations of fluctuation operators 2? m (oj) and 



Jo(r,w) [12|. The scattering intensity is evaluated 
by the first-order correlation C^(x) — (E~(x)E + (x)) 
at x — (r,uj), where r is a position outside the 
material and the scattering frequency ui differs from 
the pump frequency w m . The two-photon coinci- 
dence intensity is evaluated by the second-order corre- 
lation C ( - 2 \x u x 2 ) = (E~ (xi)E~ (x 2 )E+ (x 2 )E+ (xi)) = 

(2) (2) (2) 

Cg (xi,X2) + C N (xi,X2)- Here, C s (xi 1 x 2 ) is the in- 
tensity of correlated photons satisfying the frequency 
conservation 2wi n = uij + uj 2 , while C N (xi,x 2 ) = 
C"(ji)CW(i2) is the intensity of accidental pairs and 
is written as a product of two first-order correlations [l5| . 

We consider a CuCl film with thickness d, and the z- 
axis is perpendicular to its surface. We consider only 
the lowest relative motion of excitons and biexcitons and 
assume sinusoidal wavefunctions for their center-of-mass 
motion (indexed by n) as seen in Fig. [Ub). We assume 
their standard parabolic energy dispersions and follow 
Ref. for other excitonic parameters. Concerning biex- 
citons, its translational mass is = 2.3m ex , and the 
binding energy is A^ x = 32.2 meV Q. From Ref. 1 171 we 
consider the effective volume as \<P\ 2 = 80 nm 3 and the 
damping rate as 7b x = 13.2 fieV. The pumping light is 
a continuous plane wave perpendicular to the film, and 
has the same frequency as the giant two-photon absorp- 
tion of CuCl except in Fig. 2(b), where Ld- ln is tuned to 
resonantly excite n = 6 biexciton state. In both cases, 
hoj- m ~ hajT — Zibx/2, where Hu>t — 3.2022 eV is the 
transverse bare exciton energy. We consider the pump 
power as I — 10 /jW (except in Fig. [3Jb)) as used in 
Ref. ||| while spectral shapes of the RHPS are not mod- 
ified by changing /. The scattering angle 9 is defined 
as providing the in-plane wavenumber feii = (uj^/c) sin 6, 
and the above correlation functions are transformed into 
fc||-space as x — (9,z,lj), while they depend on z only 
for whether the scattering fields are forward or backward 
with respect to the pump propagation. We define the 
horizontal (H) and vertical (V) directions of the polar- 
ization with respect to the scattering plane. 

Fig. [2] shows the polarization-resolved spectra of for- 
ward scattering intensity (proportional to [l8j]) at 
thicknesses of (a) 7 /xm and (b) 200 nm. The nonradia- 
tive damping rate of excitons is 7 CX = 0.5 meV and the 
scattering angle is 9 = 60° . In the case of the bulk-like 
thickness [Fig. Ufa)], the four peaks, namely, Ml, Mt, 
LEP, and HEP (lower and higher energy polariton) are 
reproduced as observed in experiments yjf, 6, 8]. The LEP 
and HEP correspond to the RHPS process, while Mt 
and Ml are caused by the biexciton relaxation to the 
transverse and longitudinal excitons, respectively. The 
frequencies of LEP and HEP depend on scattering an- 
gle under the conservation of energy and wavevector in 
bulk materials. Decreasing the film thickness, the LEP 
and HEP peaks diminish due to the relax of the wavevec- 
tor conservation. Instead, multiple peaks appear in the 
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FIG. 2: Forward RHPS spectra of CuCl film with thicknesses 
of (a) d = 7 /xm and (b) d — 200 nm. V- and H-polarized 
spectra are plotted by solid and dashed lines, respectively, 
(c) Dotted lines represent the dispersion curves of upper and 
lower branch polaritons (UBP and LBP), and vertical bars in- 
dicate the resonance frequency of V-polarized exciton-photon 
coupled modes in the film with d — 200 nm. The length 
of the bars represents the sum of radiative and nonradiative 
decay rates of the coupled modes, n denotes the index of 
the original exciton state of each coupled mode. Parameters: 
7 CX = 0.5 meV, 9 = 60°. 



LEP-HEP frequency region as seen in Fig. HJb) . These 
anomalous peaks can be interpreted by the energy level 
structure of exciton-photon coupled modes in the nano 
film [§- 11 , Til , which are shown in Fig. [2fc) as vertical 
bars with length of sum of the radiative and nonradiative 
damping rates. The scattering spectra reflect the anoma- 
lous level structure of the coupled modes, because a biex- 
citon relaxes into a coupled mode by emitting a photon (a 
peak on the lower energy side) with satisfying the energy 
conservation. It is worth noting that we can significantly 
modify resonance frequencies and radiative lifetimes of 
the exciton-photon coupled modes by designing the ma- 
terial structure, i.e., material shape, size, arrangement, 
and external environments [9j, 1 1 lL 1 1 6j j . Since the anoma- 
lous level structures depend even on propagation angle 
and polarization direction, the nano-structured materials 
have much degrees of freedom to control the frequencies, 
angles, polarizations and phase difference of generated 
entangled state. This aspect is essentially different from 
the cases of bulk crystals 0, [jj Q and of the single quan- 
tum dot systems 0, 0] ■ 

For the high-power generation of the entangled pho- 
tons, the important factors are the generation efficiency 
and also the statistical accuracy, i.e., the amount of ac- 
cidental pairs. For a pumping beam with power J, the 
signal intensity S oc C s (amount of correlated pairs) is 

proportional to J 2 , while the noise intensity N oc Cjy 
(amount of accidental pairs) is proportional to J 4 , be- 
cause an accidental pair is involved with two biexcitons. 
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FIG. 3: Thickness dependence of (a) generation efficiency 
S/I 2 and (b) performance P. Dielectric constant of outside 
is Ebg, scattering frequencies are the same as that of pump- 
ing field, and scattering is forward with 8 — 0° (black, blue, 
and green). The red lines represent backward emission from 
a CuCl film with an optical cavity, whose mode frequency is 
tuned to lot- 



This implies that, by increasing the pump power /, the 
S/N ratio decreases in contrast to the increase of S 
[6j . Here, we introduce another measure termed "perfor- 
mance" P defined as signal intensity S under a certain 
S/N ratio (/ is tuned to give this ratio). This quan- 
tity does not depend on /, and represents the material 
potential for generating strong and qualified entangled- 
photon beams. While we suppose S/N = 20 as reported 
in Ref. [fj to calculate P = S 2 /20N, the spectral shape of 
the performance does not depend on a chosen S/N ratio. 

Fig. [3] shows thickness dependences of (a) generation 
efficiency S/I 2 and (b) performance P. For simplicity, we 
suppose that the two scattering fields are perpendicular 
to the film (9 = 0°) and frequencies are u>i/2 — ^in ± + . 
The black, blue, and green lines (7 CX = 0, 0.5, and 
1.0 meV, respectively) represent the forward emission 
from a CuCl film existing in a dielectric medium with 
£b g - As seen in Fig. [3]Ja) and also in Ref. [l3j, the thick- 
ness maximizing the generation efficiency is changed by 
7 CX , and it is in the order of micrometers or more for 
CuCl crystals. 

However, as seen in Fig. [3jb), the performance de- 
creases with increasing thickness for non-zero 7 GX at 
thickness range of micrometers, because the amplitudes 
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of scattering fields decrease during the propagation in the 
absorptive film (see Fig.QJa)), and the performance does 
not increase even if 7 ex = 0. The maximum performance 
in Fig. [2Kb) is the ideal quantity, and it only depends on 
measurement conditions, such as resolutions of angle and 
frequency, but not on material parameters [Tsj | . There- 
fore, the generation efficiency (generation probability for 
one pump pulse) is limited by the statistical accuracy 
(S/N ratio) when we use bulk crystals Q. However, at 
thickness of from 50 nm to 1000 nm, it is worth noting 
that nearly ideal performance values can be obtained at 
particular thicknesses even if 7 ex is non-zero. This is be- 
cause the radiative decay is dominant in such nano films 
owing to the exciton superradiance and the entan- 

gled pairs can go outside the film without decreasing the 
amplitude. Therefore, nano films generally show a high 
performance, and this rapid decay is also desired for the 
high repetition excitation, which also recovers the signal 
intensity with keeping the S/N ratio A good per- 
formance is obtained only when the resonance energy of 
exciton-photon coupled mode is just equal to a half of 
the biexciton energy. This condition appears as an oscil- 
latory behavior in this crossover thickness region, because 
the center-of-mass of biexcitons is confined and the en- 
ergy levels of coupled modes are modified by changing 
the thickness. 

However, as seen in Fig. [3](a), the generation efficiency 
of nano films is much lower than that of bulk crys- 
tals. In other words, a strong pump power is required 
to achieve a sufficient signal intensity at the thickness 
range of nanometers. This problem can be overcome by 
controlling the light-matter coupling by using an optical 
cavity. Avoiding essential modifications of the biexciton 
level scheme and considering an existing sample, we treat 
a cavity with a low quality factor (Q-factor) as reported 
in the experiment [19j, namely a CuCl film in a dis- 
tributed Bragg reflector (DBR) cavity composed by PbF2 
and PbBr2. The DBR reflectors are considered by the 
background dielectric function eh g {r, lu). The red lines in 
Fig. [3] represent the backward emission from the cavity 
structure, where the cavity mode frequency is tuned to 
w t, 7cx = 0.5 meV, and the periods of the incident- and 
transmitted-side are 4 and 16, respectively. This system 
corresponds to the weak bipolariton regime, where the 
energy splitting between polariton and biexciton levels 
are small compared to their broadening. This situation 
is in contrast with Ref . [20J, where the strong enhancement 
of entangled-photon generation from a quantum well in a 
cavity with higher Q-factor has been discussed based on 
the biexcitonic cavity-QED picture or the strong bipo- 
lariton picture. As seen in Fig. [3](a), since the biexcitons 
are effectively created, the generation efficiency is signif- 
icantly enhanced at thickness of nanometers, and it is 
larger than the maximum value in the previous situa- 
tion (7 CX = 0.5 meV, blue line). The enhancement also 
occurs when the energy of polariton (exciton-photon cou- 



pled mode) is equal to a half of biexciton energy, which 
is consistent with Ref. Comparing to the previous 
data (blue line), the period of the oscillation is doubled, 
because the RHPS process involving biexcitons with odd- 
parity center-of-mass motion is forbidden in an one-sided 
optical cavity. On the other hand, as seen in Fig.[3fb), at 
thickness of micrometers, the performance is suppressed 
compared to that by the bare CuCl film. This is be- 
cause of the multiple reflection inside the cavity, and the 
scattered fields nonradiatively decay during the propaga- 
tion. In contrast, at thickness of nanometers, especially 
at 80 nm, the performance almost keeps the ideal quan- 
tity, because of the enhancement of excitons' radiative 
decay rate by the optical cavity or the exciton superradi- 
ance 0, [l(| . In this way, for the optimum condition, we 
can highly control both the generation efficiency and the 
performance by manipulating the confinement of biexci- 
tons and the level structure of the exciton-photon coupled 
modes in a nano film embedded in an optical cavity. 

In summary, semiconductor nano-films have much de- 
grees of freedom to control generated states of entan- 
gled photon pairs. They show a high performance from 
the viewpoint of statistical accuracy, and the trade-off 
problem in bulk materials can be overcome by the en- 
hancement of radiative decay of excitons. By using a 
cavity under the weak bipolariton regime, the generation 
efficiency can be enhanced with keeping the high perfor- 
mance. We believe that our results make a breakthrough 
of entangled-photon generation with both a high gener- 
ation efficiency and the ideal statistical accuracy. 
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for JSPS Research Fellows. 



SUPPLEMENTARY MATERIAL 



Hamiltonian 



We consider the Hamiltonian on the main paper as 



H = H r 



Hr, 



(A.l) 



where H ex describes the excitonic systems, -ff ros repre- 
sents a reservoir for the nonradiative damping of exci- 
tons, H mt is the exciton-photon interaction, and H om 
describes the electromagnetic fields and a background di- 
electric medium, which is the Hamiltonian discussed by 
Suttorp et al. in Ref. |2l| and also used in Ref. |12j. We 
consider the Hamiltonian of excitonic system as 



(A.2) 
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where 6 M is the annihilation operator of an exciton in 
state /j, and w M is its eigenfrequency. We treat the exci- 
tons as pure bosons satisfying 



(A.3a) 
(A.3b) 



and their non-bosonic behavior is described by the 
exciton-exciton interaction, the second term in Eq. (IA.2I) . 
The reservoir H [os is written as 



/>oo 

.. Jo 



+ h + b l] WW^o) + g ;(n)dl(n)] }, 

(A.4) 

where d^(f2) is the annihilation operator of harmonic os- 
cillator with frequency fl interacting with excitons in 
state fi, and g^O) is the coupling coefficient. The os- 
cillators are independent with each other and satisfy the 
following commutation relations: 

[d^Q), d\, (/?')] = 6^,6(f2-n'), (A.5a) 
[d M (/?) J d„/(fl / )] = 0. (A.5b) 

Further, H mt is simply written as a product of electric 
field E(r) and excitonic polarization P cx (r): 

H int =~Jdr Pcx(r) • E(r). (A.6) 

Here, the excitonic polarization is represented as 

Pex(r)=^7' jU (r)6 M +H.c., (A.7) 

where the coefficient 'P^v) is represented by the exciton 
center-of-mass wavefunction g^ir) and unit vector e M of 
polarization direction as 



V^r) = Me^g^r) 



(A.8) 



The absolute value of M can be estimated by the 
longitudinal-transverse (LT) splitting energy Z\lt = 
|M| 2 /e £bg of excitons and by background dielectric con- 
stant £be- of excitonic medium. 



Equation of motion 

According to Ref. [12f or the QED theories of dispersive 
and absorbing media [2ll423l j , the equation of motion of 
the electric field is derived in frequency domain as 



V x V x E + (r,co) - — e h Jr,uj)E + (r,uj) 
= ifi u J (r, uj) + fi a u 2 P+(r, uj). 



(A.9) 



Here, £\, g (r,uj) is the dielectric function of the back- 
ground medium. Jq(t,lj) describes the fluctuation of 
electromagnetic fields and satisfies 



J (r,w),{J (r>'*)}t 



Jo(r,u)),J (r',-w') 



= S(u - uj')8{r - r')^-lm[e h Jr,uj)}l. (A.10) 

In the same manner as in Ref. [l2| , under the rotating- 
wave approximation (RWA), we obtain the equation of 
excitons' motion in frequency domain as 

[Hup -tiw - i7c X /2] bn{u)) 

dr V;(r)-E + (r,uj)+'D^uj) 



Acut 



in 



(A.11) 



where 7 CX is the nonradiative damping width (defined 
by {g^Q)} in Eq. (D7) of Ref. Q), and V^uj) is the 
fluctuation operator satisfying 



= <V,/*'o(u> - uj ) 



2tt 



(A.12) 



The last term on the right hand side of Eq. (|A.11I) is the 
nonlinear term due to the exciton-exciton interaction. 
Here, we define a new operator 



(A.13) 



[1,1* 



which annihilates a biexciton (two excitons) in state A or 
describes a two-exciton eigen state B\\g) by applying it 
to the matter ground state |g). The coefficient F\ ^ M is 
invariant by the exchange of two exciton indices as 



Further, it is ortho-normal as 
and also has a completeness 



(A.14) 



(A.15) 



(A.16) 



From the excitonic Hamiltonian (|A.2I) . the coefficient 
F\,n,v an d eigen frequency Q\ of biexciton eigen state 
A should satisfy 

(A.17) 
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By using Eqs. (IA.14[) and (|A.16[) . we can rewrite 
Eq. (|A~13|) as 



(A.18) 



Therefore, from this relation and Eq. (|A.17[) . we can 
rewrite Eq. (|A.11|) as 

[hUfj, - Huj - i7cx/2] b^(uj) 

dr V;(r)-E + (r,uj)+V^Lu) 



x / doj' {b u {uj' -uj)Yb x (uj'). 



(A.19) 



On the other hand, the equation of motion of B\ is de- 
rived in frequency domain as 

(hf2 x - Huj)B\(uj) 

/OO 
da/ {hujjj - &*/)&„ (u/J&^w - uj'). 
,, ., -oo 

(A.20) 



In principle, the biexciton-resonant hyperparametric 
scattering (RHPS) process is described by equations of 
motion (|A.9[) . (|A.19|) . and (|A.20I) . and commutation re- 
lations (|A.10|) and (|A.12[) . However, in the actual calcu- 
lation, we use the following approximation. 



we simultaneously solve Eq. (IA.9I) and 

[hujfj, - Huj - i7cx/2] b^(uj) 

dr V; i (r)-E + (r,uj)+V^uj) 

+ y^toju + Huj v - HQ\)F\ tll>v 

X,v 

POO 

x / duo' {bW(uj'-u,)yB x (uj'), 



(A.22) 



and then represent E + (r,uu) by the fluctuation opera- 
tors Jo(t*,w) and V^{uj). The calculation procedure is 
straightforward by u sing the Green's function technique 
(see section □ or Ref. [12[). 

For the calculation of B\(ui), we suppose that the biex- 
citon amplitude does not decrease by the RHPS process, 
because its contribution is small compared to the incident 
light. Under this approximation, by phenomenologically 
introducing a damping constant 7b x , the biexciton am- 
plitude is obtained from Eq. (|A.20[) as 



HQ\ — Huj — ■ ; 

duu' (^-^(SCD^'JXSWCw-w')), 

(A.23) 



where (bv(u)')) can be calculated from Eqs. (|A.9[) and 
(|A.21j) by considering an incident light beam as a homo- 
geneous solution of Eq. (|A.9|) . Under the weak bipolari- 
ton regime, where the coupling between polariton and 
biexciton is small enough compared to their broadening, 
the approximated expression (|A.23[) of biexciton ampli- 
tude is sufficient for the discussion of RHPS process. 



Approximation for RHPS process 

We suppose that a coherent light beam resonantly ex- 
cites biexcitons and their amplitude is large enough com- 
pared to the vacuum fluctuation. In this case, if we do 
not consider the other higher processes, the biexciton op- 
erator in the nonlinear term of Eq. (|A.19j) can be replaced 
by its amplitude B\(u>') = (B\(u)')). Further, we replace 
b w (uj' — uj) in the nonlinear term by bu{oj' — uj), which 
satisfies the linear equation 



[huy. - Huj - i 7ox /2] S«(u/) 



dr V;(r)-E+(r,uj)+'D^ 



UJ 



(A.21) 



Simultaneously solving this equation and Eq. (IA.9I) . 
bp~\uj) can be expressed by the fluctuation operators 
Jo(r,u) and T>^(uj). Under the above approximation, 



Model of biexcitons 

Although Fx ^ and f2\ should be determined from 
(|A.17|) for given nonlinear coefficient V^^-y y in princi- 
ple, we instead suppose them from the experimental re- 
sults. This treatment is useful because we already know 
many parameters of the lowest biexciton level in CuCl 
by the longstanding experimental and theoretical studies 



It is well known that the lowest level of biexciton in 
CuCl is singlet and has zero angular momentum because 
of the exchange interaction between two electrons and 
between two holes Q . Since we suppose the resonant ex- 
citation of the lowest level, we only consider the lowest 
relative motion of biexciton. Further, according to the 



RHPS experiments in Ref. |24|, the lowest biexciton state 
mainly consists of Is excitons, and the contribution from 
the higher exciton states was estimated in the order of 
1(T 4 . Therefore, we consider only Is relative motion of 
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excitons, which has a degree of freedom of polarization 
direction £ M = {x, y, z}. The coefficient Fx.^.u is propor- 
tional to the polarization selection rule as 



(A.24) 



Considering the relative motion ty(r) of two excitons in 
the lowest biexciton level, the coefficient is written as 



F x ,w = <W / dr J dr' &(r')G x (r) 9 ;(r + r')gt(r), 

(A.25) 

where g m (?") and Gi(r) are center-of-mass wavefunctions 
of excitons and biexcitons, respectively. Here, we sup- 
pose that the Bohr radius of biexciton is much smaller 
than the crystal size, and then we approximate the above 
expression as 



F x ~ S X ^ U <I> J dr Gx(r) g*(r) g* v (r), (A.26) 
where # is defined as 

<P = J dr W(r), (A.27) 



and \<P\ 2 represents the effective volume of the lowest 
biexciton state. Therefore, we can separately consider 
the relative and center-of-mass motions of biexcitons. 
\<P\ 2 has been estimated by an experiment [13], and was 
also used as a parameter in calculation 25|. 



Observables 

For the one-photon scattering intensity at r with po- 
larization direction e and frequency u> (not equal to pump 
frequency w m ), its intensity is proportional to the first- 
order correlation function 

C (1) (r,e,w) = e ■ (E~(r, u)E + (r, u)) ■ e. (A.28) 

This function can be calculated by the commutation re- 
lations (|A.10[) and (|A.12j) . The two-photon coincidence 
intensity between (rj.,ei,o;i) and (r 2 , e 2 , co 2 ) is propor- 
tional to the second-order correlation function 



C ( - 2} (r 1 ,e ll uj 1 ;r 2 ,e 2 ,uj 2 ) = (ei • E (r 1 ,uj 1 )e 2 ■ E {r 2l uj 2 )e 2 • £? + (r 2 ,w 2 )ei • £; + (ri, wi)) 

= C y s ; (ri,ei,wi;r 2 ,e 2 ,w 2 ) + C y N ' (n, e 1 , w x ; r 2 , e 2 , ui 2 ), (A. 29) 



where we neglect the interference term that appears only 
(2) 

under oj\ = uj 2 . C s represents the signal intensity or the 
number of correlated photon pairs, which satisfies the en- 
ergy conservation uj\ + uj 2 = 2uj- m . On the other hand, 

(2) 

C y N ' appears for arbitrary pair of u>i and u> 2 . This repre- 
sents the accidental coincidence of emitted photons from 
independent biexcitons, because it is just the product of 
two first-order correlation functions 



we can rewrite Eq. (|A.9[) as 



E+ (r, u) = E+ (r, lu)+^lo 2 J dr' G(r, r', oj)-P+(r', u), 

(A.32) 

where ^g~(r,o;) is the electric field in the background 
(-Hem) system, and it is defined as 



C%'(r 1 ,e 1 ,u 1 ;r 2 ,e 2 ,u 2 ) 



CW(r 1 ,e 1 ,a; 1 )C«( 



r 2 ,e 2 ,uj 2 



(A.30) 



E+(r, w) = iiiQUJ I dr' G(r, r', uj) ■ J (r', w). (A.33) 



Further, from Eq. (1A.10[) . E^(r,ui) satisfies [ll 



Solving wave equation 

Here, we explain how we simultaneously solve 
Eq. (|A.9[) and equation of exciton motion. By using the 
dyadic Green's function satisfying 

VxVxG(r,r',w) j£b g (r, w)G(r, r' , lu) — S(r — r')l, 

(A.31) 



E+(r,tj),Eo{r',w')\ = [E+ (r , w) , E+ \r' , -a/) 
= 5{u -o/)^W,r>) - G*(r,r',o;)]. (A.34) 

l27T 



The expression of G(r, r', to) in planar system (dielectric 
multi- layer) is already known [26j. 

Substituting Eq. (|A~32l) into Eq. (|A~22|) . we obtain the 
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equation set for exciton operators as 



A,// 



x / oV {8W(w' - W )}tB A (o;'), 



(A.35) 

where the coefficient in the left-hand side is defined as 
M = [H -tvjj- i7ox/2] 5 M , M < 

-Mo^ 2 / dr /" dr' V;(r)-G(r,r',uj)--P^(r'). 

(A.36) 



(A.37) 

This simultaneous linear equation set is solved by the 
inverse matrix W(w) = [S(w)] -1 , and the commutation 
relation of SjL (w) is derived in Ref. [HI as 

S«H,{$£V(cO}t 

= *( w - w') A - W^H] , (A.38a) 

(A.38b) 



S«(u,),6«(-u/) 



= 0. 



Further, Eq. (|A.35j) is rewritten as 

/>OC 

xV, / dw' {6W(c'- w )}tB A (u;'), 



(A.39) 

and then the expression of electric field is obtained by 
substituting this into Eq. (|A.32|) . 



Intensity estimation 

We have estimated the scattering intensity, signal S, 
and noise N of RHPS process from the experimental data 
in Ref. @. For a bulk crystal, the reported maximum 
count of entangled pairs (sum of HH and VV) is 50 in 
300 seconds and the width of the pulse is 2 ns. Since the 
time resolution is 300 ps Q , the number of pairs in unit 
time is roughly estimated as 

- x 50 x 2nS x — !— -0.5 s- 1 . (A.40) 
2 300 ps 300 s y ' 



The reported signal-to- noise ratio is S/N = 20, and then 
the number of uncorrelated pairs is 0.025 s _1 . The re- 
ported pump power is I — 10 /xW. For the intensity es- 
timation of our calculation results, we have suppose that 
these experimental results correspond to the our data in 
the case of CuCl film in vacuum with 7 CX = 0.5 meV and 
thickness of 7 /im, which is the optimum thickness for 
the generation efficiency S/I 2 . On the main paper, we 
consider a CuCl film with dielectric medium with £b g or 
the film with distributed Bragg reflectors. 



The last term of Eq. (|A.36[) is the renormalization due to 
the exciton-exciton interaction via the electromagnetic | 
fields. Further, the equation of operator bfi'(u)) in the 
linear regime is also rewritten as 

£ WO^M =j dr V;(r)-E+(r,uj)+t>,(u). [2 
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